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The calculation of the optical properties of doped semiconductor quantum wells is an intricate 
many-body problem because of the dynamical response of the excess carriers to the photo- 
generated valence band hole. At low densities, however, where the main effect of the dynamical 
response is the formation of trions, a simple cluster meanfield approximation can be effectively 
employed to calculate the optical susceptibility. 

I. INTRODUCTION 

The interaction of an optically generated exciton with excess carriers in the conduction band has been 
the subject of extensive experimental and theoretical investigations. Especially the high density regime 
attracted much attention and is qualitatively understood Dynamical screening and Pauli blocking 

prohibit the formation of bound eh Q states, the exciton is therefore unstable and its spectral weight 
distributed among eh scattering states. Residual correlations between the valence band (VB) hole and 
conduction band (CB) electrons are still important, however, and, depending on the CB electron to VB 
hole mass ratio, give rise to a more or less pronounced Fermi edge singularity 5], for instance, in optical 
absorption spectra. 

With the seminal experiment of Kheng et al. Q , which showed that in a weakly n-doped semiconductor 
quantum well (QW), photogenerated eh pairs give rise to trions, i.e., bound eeh states comprising two 
CB electrons and one VB hole, the center of interested shifted to low densities. As a result of numerous 
theoretical and experimental studies [7l. 1^. I9I [Tol fTTl [T^. [T^ [T^. [Tfil Il7l EM fl9l Eol I2H I2H l^. the 
importance of eeh states in the low density regime is now unambiguously established and a qualitative 
understanding of the low density regime is also starting to emerge. 

Both the Fermi edge singularity and the trions are consequences of the dynamical response of the 
CB electrons to the sudden appearance of the photogenerated VB hole. A typical lowest order process 
is depicted in Fig. £Q. The photogenerated VB hole scatters from one momentum state to another 
and simultaneously excites a virtual ee pair (recall our notation 3) to compensate for the momentum 
transfer. At low densities, the Coulomb interaction is strong enough to correlate the VB hole with two CB 
electrons (the photoexcited CB electron and the CB electron from the ee pair) and a (negatively charged) 
trion appears. At higher CB electron densities, an arbitrary number of ee pairs is excited, giving rise 
to, among others, the screening of the Coulomb interaction. At some density, the Coulomb interaction 
is then too weak to support trions and/or excitons, but the residual correlations might produce a Fermi 
edge singularity. 

The close connection between trion formation at low densities and the appearance of the Fermi edge 
singularity at high densities has been verified experimentaly |20ll2lU22| but, so far, a unified theoretical 
description of the full density dependence of the optical susceptibility, based, e.g., on the qualitative con- 
siderations of the previous paragraph, is still missing. To construct such a theory is a rather formidable 
task. Here, we focus therefore on the low density regime, where the cluster meanfield approximation 
(CMFA) |2J| can be effectively used to calculate the optical susceptibility. ^3] In contrast to early theo- 
retical studies of trion states 0, H, El El > the CMFA describes trion and exciton states simultaneously. 
Moreover, it explicitly accounts for a (low) density of excess carriers. [In this sense it is closely related 
to the density matrix method of Esser et al. an d the exciton-electron T-matrix model of Suris et 
al. ^H-] The high density regime, however, cannot be systematically reached, despite the earlier opti- 
mistic assessment of the author. 



2 




FIG. 1: The left panel illustrates a typical dynamical response of the CB electrons to the photogenerated VB hole 
in an idealized n-doped semiconductor (/j, e is the chemical potential for the CB electrons). The right panel shows 
the corresponding Feynman diagram. The dashed forward running line denotes a VB hole, whereas the forward 
and backward running solid lines depict, respectively, a CB electron and CB hole. Multiple scattering, indicated 
by the thin stripes attached to the participating lines, correlates the intermediate eeh cluster to a trion. 

II. CLUSTER MEANFIELD APPROXIMATION (CMFA) 

A. Optical susceptibility 

We consider for simplicity an idealized strictly two-dimensional quantum well with two isotropic 
parabolic bands. With an appropriate choice of single particle states, our formalism can be of course also 
applied to more realistic QW models, which include, e.g., valence band mixing and finite band off-sets. 
We assume a single hole in the VB and a low concentration n of electrons in the CB. The charge carriers 
are coupled through the Coulomb interaction and the total Hamiltonian is given by 

H = H + V 

ka kpqar 

where e c (k) = E g + k 2 /2m e and e v (k) = k 2 /2m,h are the CB and VB dispersions (h — 1), respectively, 
and V(q) — 4ire 2 /eoq is the Coulomb interaction in standard notation. 

If the momentum dependence of the dipole matrix element r 7 is ignored, linear response theory gives 
for the optical susceptibility [T stands for fci] 

xH = 2|r 7 | 2 ]T((&_ iai ;4&L 2 }} w (2) 
12 

where((6_iai; a\b*_ 2 ))u denotes the Fourier transform of the two-time correlation function 

«(&_ iai )(;); (4&L 2 )(0))) = -»e(t)([(6_ l0l )(t), (4&L 2 )(0)]>. (3) 

Here, the thermodynamic expectation value is defined by (...) = Z^ 1 Tr(e~ l3 ' H ...), with Z = Tre~ l3 ' H , 
and TL = H — fi e N e . The imaginary part of determines, apart from a constant factor, the optical 
absorption I a i, s (uo) ~ —Imx{<-o). The overall factor of two in front of the sum on the rhs of Eq. comes 
from the spin. 

B. Dyson- type equations 

The CMFA, originally developed to attack the nuclear many-body problem, rearranges the equations 
of motion for two-time correlation functions into an hierarchy of Dyson- type equations. It is thus closely 
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related to the Mori memory function method, although the mathematical details are quite different. Using 
CMFA techniques [24j, we can derive a formally exact Dyson- type equation for the eh pair propagator 
P(12;w) = ((6_iai; a\b^_ 2 ))^. Keeping in mind that the spin configurations of the various correlation 
functions are fixed, because of the spin-independence of the Coulomb scattering, we suppress the spin 
variables and write: 

[u + it] - Z(1)]P{12; u) = N(1){S 12 + ^[M s *(13) + JAf (13; w)]P(32; w), (4) 

3 

with 

M st (12) = N- 1 (l)N-\2)([[b- 1 a 1 ,V],4bU}) (5) 

8M{l2-u) = N- 1 (l)N- 1 (2)(([b. 1 a 1 ,V};[V,4bl 2 })) u (6) 

N(l) = {[b- iai ,4bl 2 ]) (7) 

Z(l) = €c (l) + e v (-l) (8) 

Applying the CMFA technique to the correlation function defining <5M(12;w), and so on, we obtain 
an hierarchy of Dyson-type equations for correlation functions with an increasing number of particles. 
Clearly this set of equations has to be truncated and it is the truncation which restricts the CMFA to 
low CB electron densities. 

We now briefly recall the truncation procedure adopted in Ref. |l2j. The vanishing VB hole concen- 
tration implies (b\bi) — 0, i.e., N(l) = 1 — (a\ai). By assumption, the CB electron density is also small, 
i.e., N(l) ~ 1. Thus, in leading order in the CB electron density, we ignore Pauli blocking. In the low 
density regime, we also neglect screening and single-particle selfenergy corrections and Eq. Q reduces 
to 

P(12; u) = P (0) (12; u) + ^ P (0) (13; ui) [M st (34) + 5M(34; w)]P(42; uj), (9) 

34 

with a bare eh pair propagator, (12; lu) = 5i2[lu + ir\ — Z(l)]" 1 , and a static part of the eh pair 
selfenergy M s *(45) = — u(4 — 5). The dynamical response of the CB electrons, i.e., the creation of virtual 
ee pairs, is encoded in the dynamical (resonating) part of the eh pair selfenergy, which, in the low density 
regime, becomes 

5M(12;w) = ^\(5>(7)(P(l,4,4 + 5,5- 1|2, 6, 6 + 7, 7 - 2; w) 

4567 

- R(l-5, 4,4-5, -1|2,6,6 + 7,7- 2;w) 

- P(l,4,4+ 5,5 -1|2- 7,6,6 -7,-2;w) 

+ P(l-5, 4,4-5, -l|2-7,6,6-7,-2;w)). (10) 

We defined an eightpoint eeeh function P(1235|5678; ui) = ((64030201; a\a\ajb\)) ^ , for which, as indicated 
above, we again derive a Dyson-type equation, symbolically written as R(uS) — R(°^(ui) + R(°\u))[K st + 
SK(lu)]R(uj). [Assuming that the spin configuration of the photogenerated eh pair is (Tl), the spin 
configuration of the eeeh cluster defining Pis (j era J.); a is then summed over in Eq. I|10|l .] In the dilute 
limit, we need <5M(12; ui) only to o(n). Using an e-line expansion for <5M(12; ui), we can show (i) that the 
term 8K(ui) gives rise to o(n 2 ) contributions to o"M(12;w) and is thus negligible at low densities, and (ii) 
that, in leading order in the CB electron density, the eightpoint eeeh function appearing in Eq. l|lUfl in 
fact factorizes into 

P(1234|5678; w) = <5 37 .fe(3)G3(12|56; uj + e c (3)), (11) 

with / e (3) =< 0303 > and a sixpoint eeh function G5(12|34;w) = ((b5^i^2Ct2ai',a 3 a\bl_ 3 _ 4 )) LJ , which 
satisfies 

[u + it] - Z 5 (12)]G 5 (12|34; ui) = 7(12|34) + ^ V 5 (12|67)G 5 (67|34; uj), (12) 

67 

with 



V 5 (12|67) - v{7 - 2)J 6 , 1+2 _ 7 - v(l - 6)S 2< 7 - v(2 - 7)<J 6 ,i, 



(13) 
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and 

Z 5 (12) = e c (l) + e c (2)+e t ,(5-l-2). (14) 

Thus, at low densities, the hierarchy of Dyson-type equations can be truncated at the level of a sixpoint 
eeh function and, in leading order in the CB electron density, the optically generated eh pair is only 
coupled to an eeh cluster. The groundstate of the eeh cluster has antiparallel CB electron spin. Focusing 
on the eeh groundstate, we fix therefore the spin configuration of the eeh cluster defining G to (Til)- 
For this spin configuration, J(12|34) = 613624, and the spin eventually leads to the overall factor of two 
in Eq. J2J) for the optical susceptibility. 

Using the spectral representation for the sixpoint eeh function, 

G5(12 |34;,) = S ^a%M , W 
' uo + irj — \l n (o) 

we finally obtain for the dynamic part of the eh pair selfenergy 

^(34;. ) =E/ e (5) JT |fg^ ) (16) 

on 

where all multiple scattering events within the eeh cluster are now encoded in a vertex function 

3 5 »(3) = «(5 - 7)[#5«(37) - *Bn(7, 3 + 5-7)], (17) 
7 

given in terms of (normalized and complete) momentum-space eeh wavefunctions v I / 5„(37), which satisfy 
an eeh momentum-space Schrodinger equation: 

[Q n (5) - Z6(12)]*Bn(12) - V 5 (12|34)* 5n (34) = 0. (18) 

34 

Here, 'n' and '5' depict, respectively, internal quantum numbers and the (center-of-mass) momentum of 
the (propagating) bound eeh cluster. 

To calculate P (12; ui) it is advantageous to split Eq. 10 into two coupled integral equations with, 
respectively, M st (= —v) and 6M as kernels: 



P st (12;uj) = P^(12;uo)-^P^(l?,;u J )v(i~A)P st (A2;uj) (19) 

34 

P(12;oj) = P st (12; to) + J2P st (13; lo)SM(M)P(42;lu) (20) 

34 

The first integral equation is readily solved using again a spectral representation 

^ u) + vi] — hj v 

where the (normalized and complete) momentum-space exciton wave functions $1/(1) obey an eh 
momentum-space Schrodinger equation (Wannier equation) 

[E v - Z(l)]$ v (l) + v(l - 2)$„(2) = 0. (22) 
2 

The solutions of this equation are analytically known. pE| 

The second integral equation cannot be solved exactly. In order to avoid heavy numerics, we employ 
for its solution a simple approximation scheme based on the "exciton representation", i.e., on a bilin- 
ear expansion in terms of the solutions of Eq. J52J. Writing P(12;ui) = „ -P Z y^(^j) < I ) l /(l)* ( I ) ^(2), we 
transform Eq. I|20|l into an infinite set of algebraic equations for the expansion coefficients, 

[u + if! - E v ]P uli (u) = 5^ + Y^ 8M vX (u)P Xlx (w), (23) 
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FIG. 2: Graphical illustration (not a Feyman diagram) of the dynamical part of the exciton selfenergy 8Mx{oj)- 
The full circles denote the exciton-trion coupling constants u£ T and [u^]*. 

where 

^M=s/.w M+t ^: n . (5) . <*> 



5n 



and = $*(1)55„(1). In the vicinity of the exciton resonance, i.e., for energies Li close to the 
exciton groundstate energy Ex, we expect the exciton and trion groundstate to be dominant. Thus, we 
solve Eq. I|23[l by employing one-pole approximations: First, we keep only the term corresponding to the 
exciton ground state v — X and neglect the coupling to excited exciton states v ^ X. Next, assuming 
that the main modification of the exciton ground state comes from the groundstate of the eeh cluster, 
i.e., from the trion denoted by n = T, we restrict the sum in Eq. (|24() to the term n = T. 

Let us emphasize, the one-pole approximations are by no means mandatory and various improvements 
are possible. For example, if the exciton binding energy is small, i.e., if the exciton groundstate is not 
well separated from the first excited state, the first excited state has to be kept. This gives rise to two 
coupled algebraic equations (two-pole approximation). It is even conceivable to treat all eh scattering 
states in a similar way, introducing an "effective" excited state, which carries the spectral weight of the 
whole continuum. Excited trion states can be also taken into account by simply extending the sum in 
Eq. (|24() to the respective states. 

The theoretical work by Esser et al. indicates that these refinements are indeed necessary to obtain 
good quantitative agreement between theory and experiment. For the purpose of demonstration, however, 
we neglect here scattering states and stick to the simple one-pole approximation described above. The 
optical susceptibility in the vicinity of the exciton resonance is then given by 

X H ~ — 2 1 SXI \ M . y (25) 
with an exciton selfenergy (see Fig. [2J) 

5M x H=^/ e (5) |M 7 |2 (26) 

uj + vr\ + e c (5) - O t (5) 

and an exciton oscillator strength sx = r 1 J^i ^x(l)- Note, the one-pole approximation is conceptually 
close to the T-matrix model considered by Suris et al. in their analysis of absorption and reflection 
spectra of modulation-doped quantum wells. [l5| 

To proceed further, we assume that the CB electron density is such that the Fermi function in Eq. 
(|26|l can be replaced by a Boltzmann function parametrized by the CB electron density n. Measuring 
energies in units of the 2D exciton binding energy 4i? (R is the 3D exciton Rydberg) and length in units 
of the 3D Bohr radius as, we define x(^) = x(4^ + Eg)/4R, with ui = [uj — E g )/4R, and rewrite the 
optical susceptibility into dimcnsionless form 

m = 'glV-* M ~ f (-V (27) 
on* uj + ir\ + 1 — dMxyu) 



Here, fj = rj/AR and 



5M x {Cj) = 5M x (4:Ru + E g )/4:R (28) 
■ Mr 
M x 



= ^(na 2 B )\uX\ 2 m, (29) 
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where (3 = 4/3 R, Mx, and Mr are, respectively, the thermal energy measured in units of AR, the exciton 
mass, and the trion mass. The exponential integral, 

r°° e-^M^y 

m = / dy — — , (30) 

7 w + it] - e T + y 

where = [f2y(0) — 2E g ]/AR is the energy of a trion at rest, i.e., the internal trion energy, which 
is the binding energy of two CB electrons to one VB hole. The dimensionless exciton-trion coupling 
constant u?f finally is given by an integral involving the wave function of the trion groundstate, the 
Coulomb interaction, and the wavefunction of the exciton groundstate. For simplicity we neglected the 
weak dependence of on the (center-of-mass) momentum of the trion. Note, for vanishing CB electron 
density n, 8Mx{&) —* and x{Cj) describes a single exciton resonance. 

To obtain the exciton-trion coupling constant u?f and the internal trion energy £t, we employ the 
variational technique originally used by Stebe and coworkers 0, to calculate various properties of 
an isolated trion. For that purpose we transform Eq. (|18(l to real space (to solve the eeh problem we 
use atomic units, which are, however, at the end of the calculation, translated into the excitonic units 
adopted above) 

[T + V}^ T (r el r e2 f h ) = £T*T(r e i?" e2 r) l ), 
with kinetic and potential energies (h = 1), 

1 -(A e i + A e2 ) - -—Afc, 



2m e 2m?! 

e\ 1 1 1 
V = ( p ^ + p ^-p =- , 

eo \r e i - r h \ \r e 2 - r h \ \r e i - r e2 \ 

separate away the center-of-mass motion, which is simply a plane wave, and make a variational Hylleraas 
Ansatz for the internal part of the trion groundstate wavefunction, 

*(s, t, u) = c lnm \lnm), (31) 

Inm 

where \lnm) = s l t m u n and u, s, and t are ellipitical coordinates for the internal motion of the eeh 
cluster. Note, since the trion groundstate is a singlet, ^(s, t, u) has to be an even function in t, i.e., m has 
to be even. The expansion coefficients are determined by a Ritz variational principle, i.e., by minimizing 
the energy functional E[c, ex] — (W|T + V r |^E')/(5 , | x I'). The variation yields an eigenvalue problem which 
is then iteratively solved. The largest eigenvalue is related to (—* &r) and the associated eigenfunction 
gives the expansion coefficients c; nm . For more details we refer the readers to Refs. H E3 and 0. With 
the Hylleraas variational function, the exciton-trion coupling constant becomes 



/ 2 1 f°° , f s , f u , 2n(s 2 -t 2 )u _^±^ r 2 1 1T/ 
- ' ds du dt—=±==L=e W[ ~]V(s,t,u), (32) 



u 



Nnk{l + a)J J J ^(u 2 -t 2 ){s 2 - u 2 ) s-t 
where k is the effective charge of a CB electron, N is the norm of 'J, and a = m e /mh 



C. Some numerical results 

Using a 2D 22-term Hylleraas wavefunction, we calculated in Ref. ^3] the optical absorption along 
the lines explained in the previous subsections. For a mass ratio a = m e /m^ = 0.146, corresponding 
to GaAs, we found u% ~ 1.0896 and e T = -1.1151. (w£ was obtained by Gaussian quadrature.) The 
binding energy (in units of AR) for one electron is therefore Wt = 0.1151 in agreement with the binding 
energies obtained by other means. 

The optical absorption in the vicinity of the exciton resonance is related to the exciton spectral function 

D x {u) = -Imx(u>), (33) 

shown in Fig. [3] for na 2 B = 0.008 and (3 — 10. The overall structure is in qualitative agreement with 
experiments, despite the simplicity of the model. Below a narrow peak at uj = —0.94, we see a broad 
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FIG. 3: The spectral function Dx(u)) for na% = 0.008 and AR/3 = 10. The effective mass ratio a = 0.146, 
corresponding to GaAs, and fj = 0.001. The left inset depicts the integrated spectral weight Cx(0j), whereas the 
right inset displays Tx(&) (dot-dashed line) and Ax(w) (solid line). The energy for which Ax(£>) crosses Co + 1 
(dotted line) defines the exciton line at Q — —0.94. From Ref. fl^l. 



absorption band with a sharp high energy edge at w = —1.1151 and a low energy tail. The narrow 
peak is the exciton resonance whereas the broad band corresponds to the absorption due to a trion. As 
shown in Ref. |l2( , the sharp high energy edge of the absorption band comes from the pole of the exciton 
selfenergy at er = —1.1151, i.e., a photon in resonance with the high energy edge creates a trion with 
momentum zero. The band has a tail at the low energy side because of the recoil of the trion, i.e., the tail 
comprises trion states with finite momentum. The right inset of Fig. [3] displays the real and imaginary 
parts of the exciton selfenergy, SMx(<j->) = Aj(w) — iTx(&)- As can be seen, trion states with small 
momenta (high energy edge) are strongly damped. The maximum of the absorption band does therefore 
not correspond to a trion with zero momentum (as in the case of vanishing CB electron concentration |llj |) 
but to a trion state with a finite momentum. The left inset depicts the integrated spectral weight up to 
an energy Q defined in terms of the cumulant Cx(&) = f_ dEDx(E). As a consequence of the one-pole 
approximation, the total spectral weight associated with the trion band and the exciton line adds up to 
one. If scattering states had been taken into account, they would also carry some spectral weight, i.e., 
in a complete theory, the combined spectral weight of trion and exciton has to be of course less then 
one. Our exploratory calculation neglected scattering states, because, at low enough CB densities, in 
particular, within the range of validity of our theory, scattering states are expected to carry almost no 
spectral weight. Indeed, the more refined calculations of Esser et al. ^3 seem to verify this assessment. 
The spectral weight, which can be associated with scattering states, is less then 5% of the spectral weight 
of the exciton groundstate. 



III. CONCLUSIONS 



We presented a detailed description of the CMFA for the optical response of a weakly n-doped idealized 
semiconductor QW. The CMFA, tailormade for low CB electron densities, reduces the calculation of the 
optical susceptibility to the solution of the Schrodinger equations for, respectively, a single eh pair and 
a single eeh cluster. Excess carriers are only treated as a reservoir for bound state formation, which 
does not, in leading order in the CB electron density, modify the properties of the eh and eeh states. 
We derived the CMFA for an idealized QW model, but it is straightforward to generalize the formalism 
to more realistic models. Furthermore, relaxing the one-pole approximations, which we solely used for 
illustration, would allow to take exciton and trion continuum states into account. At intermediate-to- 
high CB electron density, many-body medium corrections, such as Pauli blocking, screening, and single 
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particle selfenergy corrections, substantially modify eh and eeh states. Within the CMFA it is however 
not clear how to systematically account for these effects. 
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